Observation of out-of-phase bilayer plasmons in YBa 2 Cu307_,5 
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The temperature dependence of the c-axis optical conductivity o~(u) of optimally and overdoped 
YBa2Cu30^ (x=6.93 and 7) is reported in the far- (FIR) and mid-infrared (MIR) range. Below 
T c we observe a transfer of spectral weight from the FIR not only to the condensate at uj=Q, but 
also to a new peak in the MIR. This peak is naturally explained as a transverse out-of-phase bilayer 
plasmon by a model for <j(uj) which takes the layered crystal structure into account. With decreasing 
doping the plasmon shifts to lower frequencies and can be identified with the surprising and so far 
not understood FIR feature reported in underdosed bilayer cuprates. 



After many years the discussion about the charge dy- 
namics perpendicular to the Cu02 layers of the high T c 
cuprates is still very controversial. The role attributed 
to interlayer hopping ranges from negligible to being the 
very origin of high T c superconductivity Q . There is no 
agreement about the relevant excitations nor about the 
dominant scattering mechanism. The c-axis resistivity 
p c is much larger than predicted by band structure cal- 
culations. The anisotropy p c / p a b can be as large as 10 5 
and shows a strong temperature dependence, especially 
in the underdoped regime, which has been interpreted as 
an indication for non- Fermi liquid behavior and confine- 
ment Q . This strong temperature dependence is due to 
two different regimes with d,o c /dT < for T c < T < T' 
and d / o c /dT>0 for T>T', with a crossover temperature 
T' that decreases with increasing doping. There is some 
agreement as to the phenomenology that p c is described 
by a series of resistors f|-||], i.e., that different contri- 
butions have to be added, and that the sign change in 
dp c /dT is due to the different temperature dependence 
of the competing contributions. Overdoped YI^CusOz 
(YBCO) is often regarded as a remarkable exception, as 
Pd Pab is only about 50, and dp c /dT > for all T > T c . 
It is an important issue whether a sign change in dp c /dt 
at low T is really absent or only hidden by T c being 
larger than a possible T', i.e., whether overdoped YBCO 
follows anisotropic three dimensional (3D) or rather 2D 
behavior. 

The c-axis optical conductivity <7i (u>) of YBCO shows 
several remarkable features |^||: (1) It's very low value 
compared to band structure calculations, reflecting the 
large p c . (2) A suppression of spectral weight at low 
frequencies already above T c in underdoped samples re- 
ferred to as the opening of a 'pseudogap' (which agrees 
with the upturn in p c ). (3) The appearance of an intrigu- 
ing broad 'bump' in the FIR at low T in underdoped 
samples. (4) In overdoped YBCO, the spectral weight 
of the superconducting condensate is overestimated from 
cri(u>) as compared to microwave techniques ]To| . 

In this letter we suggest that most of the above men- 
tioned issues can be clarified by modelling the cuprates 
or in particular YBCO as a stack of coupled Cu02 layers 



with alternating weaker and stronger links. This multi- 
layer model fits the measured data at all doping levels and 
at all temperatures. A similar model was proposed for the 
superconducting state by van der Marel and Tsvetkov 
pi]. A transverse optical plasmon was predicted. This 



model has been verified in SmLa .8Sr .2CuO4_5 |12|. We 
report the observation of this mode in the infrared spec- 
trum of optimally and overdoped YBCO and propose a 
common origin with the above mentioned 'bump' in un- 
derdoped YBCO. 

Single crystals of YBa2Cu30 a; were grown using the re- 
cently developed BaZr03 crucibles jlj| , which in contrast 
to other container materials do not pollute the result- 
ing crystals. Crystals grown using this technique exhibit 
therefore a superior purity (> 99.995 at. %) M. The 
samples had typical dimensions of 2x0. 5-0. 7 rnn? in the 
ac-plane. The O concentration was fixed by annealing ac- 
cording to the calibration of Lindemer . An O content 
of x~7 was obtained by annealing for 400 h at 300°C in 
100 bar of high purity oxygen. Annealing in flowing oxy- 
gen at 517°C for 260 h produced 1=6.93. Measurements 
of the ac-susceptibility indicate T c =91 K for x=6.93 and 
87 K for x—7. The width of the transitions were 0.2 K 
and 1 K, respectively. Polarized reflection measurements 
were carried out on a Fourier transform spectrometer be- 
tween 50 and 3000 cm -1 for temperatures between 4 and 
300 K. As a reference we used an in-situ evaporated Au 
film. Above 2000 cm" 1 the spectra are almost T inde- 
pendent. The optical conductivity cr(u>) was calculated 
via a Kramers-Kronig analysis. 

The measured c-axis reflectivity and <Ji(u>) derived 
from it are plotted in Fig. 1 for 4 and 100K (solid and 
dashed black lines). Disregarding the phonons, ai(u>) 
shows an almost constant value of about 200 f2 _1 cm _1 . 
A Drude-like upturn is only observed at low frequencies 
in the overdoped case x—7. Below T c a sharp reflec- 
tivity edge develops at about 300 cm -1 (inset of top 
panel), which had been identified as a Josephson plas- 
mon, a collective mode in a stack of Josephson coupled 
2D superconducting layers. The gradual suppression of 
<ti(uj) below about 700 cm -1 can be attributed to the 
opening of the superconducting gap. The finiteness of 
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<7i(tL>) at all frequencies reflects the d- wave symmetry of 
the gap. The increase of <j\{uj) between 700 and 1300 - 
1500 cm" 1 from 100 to 4 K comes as a surprise. The su- 
perconducting phase transition obeys case II coherence 
factors for electromagnetic absorption jig] , i.e., only a 
suppression of o~i(lu) is expected for frequencies not too 
close to 0. The difference of spectral weight above and 
below T c defined as (for T<T C ): 

w|(T,w)=8/ [<ri(100K,u/) -oi(T,u/)]dct/ (1) 
Jo+ 

is expected to rise monotonically with increasing fre- 
quency to a constant value for frequencies much larger 
than the gap. It is common practice to determine the 
spectral weight of the superconducting condensate from 
this constant value. However, in YBCO6.99 Homes et al. 
fiof reported a discrepancy between u>a determined from 
either this optical sum rule (2050 ± 150 cm" 1 ) or the mi- 
crowave surface reactance (1450 ±50 cm" 1 ). To account 
for this difference the existence of a very narrow normal 
carrier Drude peak with a width smaller than the low- 
est measured frequency was concluded, which contradicts 
again the microwave measurements showing a very small 
o~\(uj) |L7[ . Our data clearly indicate a non-monotonic 
behavior of loa{lo) (insets in Fig. 1, see also Ref. 0) and 
a spectral weight transfer from low frequencies to a new 
peak above the phonons. This can naturally be explained 
by the following model for o~{u>) which takes into account 
the layered structure of the cuprates. 

We devide the unit cell of YBCO into the intra- and 
inter-bilayer subcells A and B. Let us imagine, that a 
time dependent current is induced along the c-direction, 
the time derivative of which is (dJ c /dt). We define 
(dVj/dt) as the time derivative of the voltage between 
two neighboring Cu02 layers, i.e., across subcell j. Our 
multilayer model corresponds to the approximation, that 
the ratio (dVj/dt)/(dJ c /dt) is provided by a local lin- 
ear response function pj corresponding to the complex 
impedance which depends only on the voltage variations 
on the neighboring Cu02 layers, and not on the volt- 
ages on layers further away. Microscopically this corre- 
sponds to the condition, that in the normal state the 
mean free path along c must be shorter than the dis- 
tance between the layers, lj. In the superconducting 
state this should be supplemented with the same con- 
dition for the coherence length along c. In this sense, 
the multilayer model reflects the confinement of carri- 
ers in the 2D Cu02 layers. Let us treat the current as 
the parameter controlled by applying an external field. 
Since the current between the layers is now uniform and 
is independent of the subcell index j, the electric field 
average over the unit cell is a linear superposition of the 
voltages over all subcells within the unit cell. This effec- 
tively corresponds to putting the complex impedances 
pj of subcells in series, p(u>) = xaPa{u) + xbPb(u), 
where the Xj = lj/l c are the relative volume fractions 
of the two subcells, Ia + Ib = lc, and Pj(u>) are the lo- 
cal impedance functions within subcells A and B. This 



sum for p(lu) = [a(uj)+uj/ATri]~ 1 is very different from the 
case of a homogeneous medium, where different contribu- 
tions are additive in o-(oj) = Eoj (lu), which corresponds to 
putting the various conducting channels of the medium 
in parallel. To illustrate this, let us adopt the Drude 
model for the complex interlayer impedance. In parallel 
conduction the sum of e.g. two Drude peaks yields 



p{u) 



u p,A 



J p,B 



(2) 



where oj p .j denotes the plasma frequency, and jj labels 
the damping. This results in a single plasma resonance at 
a frequency ujp = lu^ a + uj^ b , i.e., only one longitudinal 
mode (the zero) survives which is shifted with respect to 
the zeros of the individual components. The transverse 
mode (the pole at uj — 0) is identical. Putting two Drude 
oscillators in series in the multilayer model, i.e., using 
T,XjPj has a surprising consequence. 



x B 



"p,A 



U P,B 



(3) 



Now both longitudinal modes (poles of pj) are unaffected, 
and in between a new transverse mode arises. This trans- 
verse optical plasmon can be regarded as an out-of-phase 
oscillation of the two individual components. This mode 
has been predicted in Ref. [|llj for the case of a multi- 
layer of Josephson coupled 2D superconducting layers. 
The existence of two longitudinal modes was confirmed 
experimentally in SmLao.sSro^CuO^tS jliq ]. Note that 
superconductivity is not a necessary ingredient, the op- 
tical plasmon appears regardless of the damping of the 
individual components. 

In order to apply the model to the measured reflec- 
tivity data we have to include the phonons, for which a 
separation into subcells is not generally justfied, e.g. for 
the c-axis bending mode of the planar O ions, located 
on the border between subcells A and B. Therefore we 
adopt the following model impedance 



o-j + o-ph + &M + w/iwi ' 



je{A,B} (4) 



where xa = 0.28, and xb — 1 — xa for YBCO. Note that 
this model reduces to the conventional expression for a 
homogeneous medium commonly used for high T c su- 
perconductors if we either set xa = or a a — o~b ■ The 
o~a,b{u) contain the purely electronic contributions with 
eigenfrequency c<jo = within each subcell. 



lul 



je{A,B} 



(5) 



where ui s j and ui n j label the plasma frequencies of super- 
conducting and normal carriers, respectively. All other 
contributions (phonons, MIR oscillators, etc.) are as- 
sumed to be identical in the two subcells and are included 
in a sum of Lorentz oscillators. 
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[<Jph + CTAf ] = V — (6) 

where wo,j denotes the j-th peak frequency. 

The agreement between the measured reflectivity data 
and fits using this model is very good at all tempera- 
tures (thick gray lines in Fig. 1). The strong MIR peak 
of the optical plasmon caused by the out-of-phase oscil- 
lation of the superconducting carriers in the two sub- 
cells is very well reproduced. Note that in a conventional 
Lorentz model the optical plasmon would have to be fit 
with three parameters too, u p and 7. Also our model has 
three new parameters, namely the two sets of oj s , lu u and 
7 of Eq. H for the two subcells as compared to the sin- 
gle set used within a conventional two-fluid fit. In the 
case of x=6.93 at 4 K we have uj n ,A — ^n,B = 0, leaving 
only one new parameter oj s . In Fig. 2 we plot the real 
part of the dynamical resistivity p(uj). The thick gray 
line was obtained from the full fit parameters and agrees 
with the Kramers-Kronig result. The solid line depicts 
the electronic contribution p e (w), which was obtained by 
leaving out the phonon part a p h(to) from the fit parame- 
ters in Eq. || In the multilayer model p e (u) is the sum of 
the subcell contributions Xjp e j = Xj/{o~j + o~m + w/47ri) 
(j 6 {A, B}, dashed lines), which shows that the two 
peaks in /O e (w) can be attributed to the plasmon peaks in 
the two subcells. Contrary to the conventional model, the 
different contributions are not strictly additive in <ti(cj) 
due to the inverse summation in Eq. || Nevertheless we 
can calculate an estimate of the electronic contribution 
a e (uj) from the fit parameters in the same way as done for 
p e . An estimate of only the normal electronic contribu- 
tion a en (u>) is obtained by leaving out the London terms 
oc uj^j together with a p h- The contribution arising from 
the presence of superconducting carriers is then defined 
as a es {Lo)=(j e (Lo)-(j en (uj) (see Fig. 1). 

With decreasing doping level the absolute value of 
<7i(w) decreases and therefore the optical plasmon peak 
becomes sharper. At the same time, all plasma fre- 
quencies and hence also the optical plasma mode shift 
to lower frequencies. This scenario explains the strong 
FIR 'bump' reported in underdoped YBCO Simi- 
lar bumps have been observed in other bilayer cuprates 
jil|,[l9]], but never in a single layer material. This bump 
has hindered an unambiguous separation of electronic 
and phononic contributions to <j\(oj). In Fig. 3 we show 
reflectivity spectra of underdoped samples of YBCO 
taken from Refs. (||| together with fits using the mul- 
tilayer model. Again good agreement with the model 
is obtained. The strong phonon asymmetries present in 
the underdoped samples called for a fine tuning of the 
model: the two apical O stretching phonon modes at 
about 600 cm -1 were described by local oscillators in the 
inter-bilayer subcell B, i.e., they moved in Eq. |J from 
o p h(u) to <tb(oj). The figure demonstrates that this re- 
produces the asymmetry of the experimental phonon line 
shape well, although a Lorentz oscillator was used. Sim- 
ilar fine tuning has only a minor effect on the quality of 



the fit for the data presented in Fig. 1. Comparing the 
various doping levels shows that both the bending (350 
cm -1 ) and the stretching (600 cm -1 ) phonon modes show 
strong assymetries whenever they overlap with the trans- 
verse plasma mode, but that both modes are symmetric 
if the transverse plasmon is far enough away, as e.g. in the 
case of x=7. Previously it was argued that the phonon 
spectral weight is only conserved for different T if the 
bump is interpreted as a phonon ||. However, a sum 
rule exists only for the total ai(uj), not for the phonon 
part separately. Moreover, in this scenario the width of 
the bump, it's temperature and doping dependence and 
the phonon asymmetries remained unexplained. 

Both the low frequency Josephson plasmon and the 
bump are suppressed simultaneously by Zn substitution 
which supports our assignment that both peaks are 
plasma modes. An increase of spectral weight of the 
bump with decreasing T was reported to start far above 
T c HU, but a distinct peak is only observed below T c . 
We obtained good fits for all T (not shown). As men- 
tioned above, superconductivity is not a necessary ingre- 
dient of the multilayer model, an out-of-phase motion of 
normal carriers will give rise to a peak at finite frequen- 
cies, too. Upon cooling below T c , the reduction of the 
underlying electronic conductivity due to the opening of 
a gap and the reduced damping produce a distinct peak. 

Our results imply that even the c-axis transport be- 
tween the two layers of a bilayer is incoherent, which 
agrees with the absence of a bilayer bonding-antibonding 
(BA) transition in our spectra. Using photo electron 
spectroscopy Q a BA splitting of about 3000 cm -1 was 
reported. The anomalous broad photoemission lineshape 
may explain the absence thereof in the optical data. 

In conclusion, we observed the out-of-phase bilayer 
plasmon predicted by the multilayer model. The good 
agreement of the optical data with the multilayer model 
at all temperatures and doping levels shows that YBCO 
can be modelled by local electrodynamics along the c- 
axis in both the normal and the superconducting state. 
This applies even to overdoped YBCO, one of the least 
anisotropic cuprates. Our results strongly point towards 
a non-Fermi liquid picture and confinement of carriers to 
single Cu02 layers. 
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FIG. 1. The c-axis reflectivity 7? and cri(w) above (dashed lines) and below T c (solid lines). The thick gray lines depict fits 
of R using the multilayer model and, in <7i(w), the normal carrier electronic contribution a en (u>) derived from it. The filled 
areas show a es (uj) as defined in the text. Insets: wa(4K,o>) as defined in Eq. 1 and detailed T-dependence of R. 
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FIG. 2. Real part of the dynamical resistivity p(ui) as obtained from the multilayer model (thick gray line). The solid 
line depicts the electronic contribution only, and the dashed lines the subcell contributions p e A(uj) and p e s(cj) to it. Insets: 
T-dependence of p e A(w) and p e s(w)- 
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FIG. 3. Reflectivity data as taken from Refs. [5,6] and fits using the multilayer model. In all panels the thick gray lines show 
the fit result. The mid panels show <7i(w) and the different electronic contributions to it (a es (ui): filled area, a en (u))\ dashed 
line) . The solid lines in the bottom panels show the electronic contributions to p(ui) and the separation into subcells A and B 
(dashed lines). 
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